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Abstract 

In this paper we explore some of the features of a heavy quark-antiquark pair at finite 
temperature using a five-dimensional framework nowadays known as AdS/QCD. We shall 
show that the resulting behavior is consistent with our qualitative expectations of thermal 
gauge theory. Some of the results are in good agreement with the lattice data that provides 
additional evidence for the validity of the proposed model. 
PACS: 12.38.Lg, 12.39.Pn, 12.90.+b 
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1 Introduction 

This is our second paper of a series devoted to the thermal properties of a pure gauge theory 
within a five-dimensional framework nowadays known as AdS/QCD. We explore further the 
model proposed in the first paper [I]. It is based on the following Euclidean background metric 

ds 2 = R 2 ^ (jdt 2 + dx 2 + jdz 2 ^j , h(z) = e^ z2 , f(z) = 1 - (^) 4 , (1.1) 

where t is a periodic variable of period irz T such that z T = with T, the temperature. At zero 
temperature, we have in fact the slightly deformed AdSs metric. Such a deformation is notable 
because it results in a Regge-like spectrum. This fact allows one to fix the value of c from the p 
meson trajectory. For our purposes, we use the estimate of [2] 

c«0.9GeV 2 . (1.2) 
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It is worth noting that the metric (|l.ip does not contain any free fit parameter. Moreover, 
the only dimensionful parameter of our model is the Regge slope of meson trajectories. It is 
quite unusual for QCD, where it is A QCD . Thus, evaluations of thermodynamic quantities and 
the Polyakov loop we are going to undertake can be considered as a further consistency check of 
our model. 

For the case of interest, let us briefly point out a couple of facts, 
(i) The free energy of a heavy (static) quark-antiquark pair at finite temperature is expressed in 
terms of a correlator of two Polyakov loops [3] 

F(r,T) = -Tln(L(xi)Lt(£ 2 )) +Tc(T), (1.3) 



with r = \xi — X2\ and x% being a point in Br. In (11. 3p the free energy is defined up to a 
normalization constant c(T) which is related to the infinite self-energy of the quarks. 

There is a subtle point here. In the literature F is often called the heavy quarkpotential at 
finite temperature. Apparently, such a definition discards the entropy contribution^ 

An order parameter for the confinement-deconfinement phase transition is the expectation 
value of the Polyakov loop. After the normalization of (|1.3I) . it is then 



L = exp{-^F(r = oo,T)}. (1.4) 

(ii) In discussing a Wilson line within AdS/CFT (QCD) [5], one first chooses a contour C on 
a four-manifold which is the boundary of a five-dimensional manifold. Next, one has to study 
fundamental strings on this manifold such that the string world-sheet has C for its boundary. 
The expectation value of the loop is schematically given by the world-sheet path integral 

(W(C)) = J DXe~ Sw , (1.5) 

where X denotes a set of world-sheet fields. S w is a world-sheet action. 

In principle, the integral can be evaluated semiclassically in terms of minimal surfaces that 
obey the boundary conditions. The result is written as 

(W(C)) =Y,Wne~ Sn , (1.6) 

n 

where S n means a regularized minimal area whose relative weight is w n H 

The paper is organized as follows. In section 2, we develop the framework we will work. In 
section 3, we present, on the basis of AdS/QCD, a few results on a heavy quark-antiquark pair in 
a thermal medium. In particular, we exhibit the free energy, the string tension, and the entropy 
at low temperatures as well as the Polyakov loop expectation value. We conclude in section 4 
with a discussion of some open problems. 



1 Fbr a discussion of this issue, see, e.g., [4] and references therein. 

2 The point is that the areas are divergent but the divergences are proportional to the circumference of C. 
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2 General Formalism 



Our basic approach will be as follows. It is believed that the correlator of the two Polyakov loops 
has a path integral representation like that of the Wilson loopH Given the background metric, 
we can attempt to evaluate the values of the regularized areas S n . If we then discard quantum 
fluctuations of strings, in the case of one dominant exponent we will get that the free energy is 
simply proportional to 

F = TS*. (2.1) 

Note that the regularized areas are defined up to normalization constants. Since these constants 
are in fact due to the infinite self-energy of the quark sources, we omit c(T). Moreover, w* can be 
absorbed into a proper normalization constant too. In this paper we will use the approximation 

But before going on, let us shortly pause here to gain some intuition about the problem at 
hand. As in pQ, we introduce the notion of an effective string tension depending on the fifth 
coordinate z. It is given b>@ 

<y{z) = - 2 ^](z). (2.2) 

z A 

Now consider the behavior of a string bit in the effective potential V = c(z). An important 
observation is that the form of V is temperature dependent. Indeed, a short algebra shows that 
there is a special value of temperature T x = i /-!= ~ 130 MeVjl Below T\ the effective potential 

has local extrema at 



^= sin(± arcsin 2£j , z max = z T y -Z= sin(§ - | arcsin , (2.3) 
while above this temperature the potential is just a decreasing function of z as shown in Fig.l. 
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Figure 1: Schematic representation of the effective potential below and above T±. 

3 Note that the corresponding world-sheet has now two boundaries. 
4 It follows from Eq.JZlJ. 
5 We use (j 1 . 2 p for all estimates. 
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Such a behavior clearly has some of the suspected properties of gauge theory at finite tem- 
peratureH For temperatures below T\ the string ended on the heavy quark-antiquark pair set 
at z = can not get deeper than z min in z direction because a repulsive force prevents it from 
doing so. This gives rise somewhat of a wall located at z = z min . The large distance physics of 
the string is determined by this wall. Since the value of the effective potential (string tension) at 
its minimum is not vanishing, the quark-antiquark free energy has a dominant linear term whose 
coefficient is proportional to a(z min ). So, this can indeed be interpreted as the low temperature 
phase. On the other hand, for temperatures above T\ the string can get deeper and finally reach 
the horizon z = z T . The large distance physics is now determined by the near horizon geometry. 
The crucial point is that the effective potential vanishes on the horizon. As a result, there is no 
linear term in the quark-antiquark free energy. This can be interpreted as thedeconfinedphaseH 

Now let us go back to the approach and choose a pair of the contours (Polyakov loops) living 
on the boundary (z = 0) of our five dimensional space. We set 

xi = (-§,0,0), x 2 = (§,0,0). (2.4) 

Next we want to look for static configurations of the world-sheet action. To this end, we 
make use of the Nambu-Goto action equipped with the background metric (jl.ip 

S = 2^7 / d ^\/ det • (2-5) 

There are basically two types of configurations to be considered. One type describes connected 
surfaces whose boundaries are the Polyakov loops, while another describes disconnected surfaces. 

Our first goal will be to analyze connected surfaces. In this case the world-sheet coordinates 
£'s can be chosen as £i = t and £2 = %■ With such a choice, the action takes the form 

r 

Sw= toT /r^^vT + W, (2-6) 

7T J- 2 Z 



where q = A prime denotes a derivative with respect to x. 
It is easy to find the equation of motion for z 



zz 



z" + (/ + (z'f) (2 - zd z \nh)-(\f + {z'f) zd z In / = (2.7) 
as well as its first integral 

hf 



(2.8) 

On symmetry grounds, we have z'\ x= o = 0. This allows us to express the integration constant 
c via the value of z at x = 0. So, we get 



c = a\ z=zo , (2.9) 

6 The notion of the effective potential turned out to be useful in studying symmetry breaking (phase transition) 
within field theories at finite temperature [6] . Like in field theory, our model also has the effective potential whose 
form is temperature dependent. 

7 Although the above arguments indicate that the phase transition occurs, we should caution the reader that it 
is a qualitative way of thinking about the problem at hand. 
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where a is defined by (|2.2p and z = z\ x= o. 

Next we perform the integral over [—§,§] of dx. By virtue of (|2.8h . it is given by 

P . 2 /*rm i _ 1 v i , , 2 ,o) 



v7 

where C is a contour in 2 plane. 

We look for solutions that obey the following condition z = maxz. The reason for this is 
that, on general grounds, the string ended on the quarks set at x = ±r/2 reaches the deepest 
point in z direction at x = 0. 

For temperatures below T±, there are two possibilities for C and, as a result, we have 



r 



a ) =vj° ^((~Y -lj 2 , with 0<z o <z mia , (2.11) 

f Z ° dz ( /<7\ 2 \~ 5 

/ M-J -lj , with z m ^<z <z T , (2.12) 



r 



but otherwise r is complex. After a short inspection we find that A 1 ^ is a continuously growing 
function of z on the interval [0, z min ]. Moreover, it equals to zero at z = and goes to infinity as 
Za —> z min . The function A 2 ^ is, unlike decreasing. It goes from its maximum (finite) value 
r m2x a t z o = -2 m ax to zero at z = z T (on the horizon). Thus, the second solution contributes only 
at distances smaller than rf^-L To complete the picture, we present the plots of in Fig. 2. 



r (fin) 




Figure 2: Typical graphs of rW below T x . Here T = 0.05 GeV. 

When temperature is increased, the interval [z min , 2: max ] becomes smaller and finally disappears 
at T = T\ . For T > T\ , a simple analysis leads to the picture that differs noticeably from that of 
Fig. 2. The point is that both the solutions now contribute only at distances smaller than some 
finite r max . This is illustrated in Fig. 3. 

Now we move on to the second type that describes disconnected surfaces. In this case a surface 
contains two pieces each of which has a topology of a cylinder. The cylinders are stretched from 
the Polyakov loops on the boundary to the horizon. If we use £i = t and £2 = z a s the world-sheet 
coordinates, the Nambu-Goto action is then 



S n 



2-kT 



dz^Jl + f(x) 2 , 



(2.13) 



r(fin) 
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Figure 3: Typical graphs of r^> above T\. r^> is growing from to r max , while is decreasing from 
r max to 0. We set T = 0.5 GeV. 



where a dot stands for a derivative with respect to z. The equation of motion for x is 



d_ 

dz 



0. 



(2.14) 



It is obvious that it has a trivial solution x = const that represents a straight string stretched 
between the boundary and the horizon. We will call the choice ()2.4p the solution Since 
this solution makes the dominant contribution, as seen from the integrand in (12. 13ft . we will not 
dwell on other solutions here. 

Having discussed the solutions, we can now see what happens with the corresponding areas. 
We begin with the solution n 1 ). To this end, we use the first integral to reduce the integral over 
x to that over z in ()2.6f) . Since the integral is divergent at z = due to the factor z~ 2 in the 
background metric, in the process we regularize it by imposing a cutoff e. At the end of the day, 
we have 



Subtracting the - term we find a finite result 



Si = — 



irTz, 



+ 



ttT 







20 dz 



h[ 1 



+ 



c(T) 



(2.15) 



(2.16) 



where c(T) stands for a normalization constant. 

For the second solution A 2 \ it is a little bit tricky because the minimal surface is built by 
sewing together two pieces. The first piece comes from A 00 ^ defined on the interval [0,z max ]. At 
z = £ max , it is sewn with the second piece coming from A 2 \ The integral (|2.13p is divergent at 
z = 0, so we regularize it by imposing the cutoff e as before. After subtracting the divergency, 
we get 



So 



7tTz„ 



dz 



dz 



C\ 2 

a 



(2.17) 



Note that both the areas are regularized in the same way, so the corresponding normalization 
constants coincide. 
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For the third solution, the minimal area can be read off from Eq, (|2.17p . One drops the third 
term as coming from the solution and replaces z max with z T in the remaining terms. The 
area is then 

Now let us look at the Si 's as functions of r. We begin with sufficiently low temperature. 
Since the solution A 2 ^ is defined only for distances smaller than r£^ x > it makes sense to first probe 
asymptotic behaviors near the point r = 0. A short inspection shows that Si is not bounded from 
below, while the others are bounded. This provides a bit of evidence in favor of dominance of S\ 
at small distances. Further numerical calculations show that this is indeed the case. Moreover, 
the first solution turns out to be dominant at physical distances tooll We present the plot of the 
regularized areas in Fig.4§ 




Figure 4: Typical graphs of Si. Here T = 0.1 GeV and c(T) = 0. 



From this figure it is also clear that there exists a critical distance r c such that for larger 
distances Soq becomes dominant. There is an apparent reason for this. At very large distances 
the quark and the anti-quark decouple from each other. Usually, instability occurs in models 
with dynamical quarks, where a string breaks. In the case of interest this occurs due to emission 
of closed strings (glueballs). As known, in AdS/QCD a free particle is described by a straight 
string stretched between the boundary and the horizon that is nothing but our solution A 00 ). 

When temperature is increased, two effects are seen: (1) The minimal difference between the 
functions S\ and 52 becomes smaller. (2) The critical distance r c is decreasing. For temperatures 
close to T = T\ , Si is no longer dominant at physical distances and it is time to account for the 
other solutions in the series (|1.6p . 



3 Applications 

In this section we will consider some applications of the developed formalism. We will focus on 
the cases where the approximation (|2.ip might be applicable. 

8 We mean the interval 0.2 fm < r < 2fm that is of primary importance for phenomenology. 
9 The overall constant can be fixed from the slope of the heavy quark potential at zero temperature. We use 
the estimate of [8] g ~ 0.94, here and below. 
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3.1 Free Energy at Low Temperature 

We begin with temperatures which are sufficiently smaller than T\. In this case the free energy 
can be evaluated by using (|2.ip with S* = Si 

First, we need to fix the normalization constant c(T). In doing so, we follow [7J and look for 
the small r expansion of the free energy. Since small distances correspond to small deviations 
in z direction, we need to study the expressions (|2.1ip and (|2.16p near z = 0. The asymptotic 
behavior of r(z ) is given b}0 

r = V^(l-fAJ + 0( 2l 5 ), (3.1) 

3 

where p = r 2 (i)/(2vr)2 . 

In a similar way we find the behavior of F. It is 

9 1 +c (T) + ^(37Tp 2 -l)z + O(zl). (3.2) 



2-7T/9 z 8irp 
Combining this with (|3.ip . we get 



F = + c(T) + a r + 0(r 3 ) , (3.3) 
r 

where kq = -§^p~ 2 and o"o = f c/? 2 . 

Having derived the small distance expansion of the free energy, we can now get that of the 
singlet free energ30 

F x = F-Tln9 = -— + c(T) -Tln9 + 0(r) (3.4) 
r 

and compare it with the small distance expansion of the heavy-quark potential at zero temper- 
ature. To leading order, we have the same Coulomb term as in [8]. To get an agreement at 
next-to-leading order, we choose 

c(T)=Tln9 + C. (3.5) 

Finally, we fix the value of C by matching it to the constant term of the Cornell potential [9]. 
This yields 

C = -0.25GeV. (3.6) 

Actually, our result (|3.4p shows that at sufficiently small distances the singlet free energy 
of the pair is temperature independent. So, the agreement with the lattice data [H [7] is very 
satisfactory at this point. 

Our next goal is to analyze the long distance behavior of F. As noted in section 2, large r 
corresponds to z ~ z min . So, we need to study the behavior of (12.111) and (I2.16P near z = z min . 
In this case a crucial observation is that the integrals are dominated by the upper limits, where 



10 In the remaining part of this section we omit the index specifying the solution each time when the meaning is 
clear from the context. 

n To this order, the calculation is identical to that of [8]. 
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Here we assume a pure SU (3) gauge theory. 
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they take the form J dv/y/a{\ — v ) +6(1 — v) 2 . Such an integral may be found in [TO]. At the 
end of the day, we have 

r = -whx(z min -z ) + O(l), F = —a T w In (z min — z ) + 0(1) , (3.7) 

where w = 2y / a/(fa")\ . and <r T = Aa\ . . 

This means that at long distances the free energy of the pair shows the desired confining 
behavior 

F = a T r + 0(l), (3.8) 

with the string tension 

e i-i / T 4x | T 2 

cr T = cr - II t jt4~ J ' t = 3^ sm^arcsm^J . (3.9) 

\ £7 s / 

Here cr denotes a tension at zero temperature. Explicitly, it is given by a = f^c [Sj. We have 
also introduced the critical temperature T tTs = \\f \ ~ 210 MeV obtained from the spatial string 
tension in pp. 

As expected, the string tension a T is a decreasing function of T. In Fig. 5 we have plotted 
a T /a against T/T\. 



cr T /<r 




Figure 5: String tension in units of a versus temperature in units of T\. 

We conclude this subsection by making a few remarks. 

(i) It is quite interesting that the tension shows very little dependence on temperature up to 
T m 0.8Ti ps 100 MeV. 

(ii) Unlike the coefficient a T of the linear term in the large distance expansion of F, the coefficient 
a of the linear term in the small distance expansion turns out to be independent of temperature. 

(iii) The free energy of the pair is written in parametric form given by Eqs. (]2.1ip and (]2.16l) . 
Since we do not know how to eliminate the parameter z and find F as a function of r, we present 
the result of numerical calculations. In Fig. 6 the free energy F versus r for a temperature below 
T\ is shown. 
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Figure 6: Free energy versus r at T = 0.1 GeV. 



3.2 Entropy and Energy at Low Temperature 

We have already mentioned that in the literature the entropy contributions are sometimes ig- 
nored It is therefore of great importance to address this issue. 

Having understood the behavior of the free energy, we can easily investigate the properties 
of the entropy. Since the free energy shows no temperature dependence at small distances, it 
makes sense to focus on long distances, where its behavior is given by (|3.8p . In this case a 
straightforward calculation leads to the following expression 

dF\ 1 fdlnf\ 

w) r --** r V«r)..- <310) 

Here we treat / as a function of two variables. Then, using the relations (jl.ll) and (|2.3I) . we can 
compute the entropy density of the pair 



_ sin 2 1 ^ arcsm =r 
3 — 4 sin ( 3 arcsin ^ 



(3.11) 



A closer look at this expression shows that the entropy density grows as the temperature 
increases. This is the expected result. To complete the picture, we plot the entropy density 
against T/T\ in Fig. 7. 




0.2 0.4 0.6 o.s 1 
Figure 7: Entropy density versus temperature in units of T\ . 



This is the case for AdS/CFT(QCD) too. See, e.g., 0. 
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Prom this figure, we see that the entropy density is close to zero up to temperatures of 
order 0.4 Ti ~ 50MeV. Thus, in this temperature range, the entropy contributions are indeed 
negligible. It is obvious that the situation changes drastically with the temperature growth. 
Therefore, a natural question to ask is whether the entropy contributions play a major role at 
finite temperature. To answer this question, we consider the internal energy of the pair. As 
usual, it is given by £ = F + TS. It is clear from above that £ is independent of temperature at 
small distances, where it coincides with the heavy-quark potential at zero temperature. At large 
distances, the internal energy has a complex dependence on temperature together with a linear 
growth with r. Explicitly, it is given by 



£ = £ T r, (3.12) 

with the tension 



2/1 a ^a\^ T 




o sin 5 arcsm ■ 2 

"t|1 + K- ^jvl • (3.J3) 

3 — 4 sin (3 arcsin ^ 



The difference with the string tension (|3.9p is due to the second factor that is nothing but 
the entropy contribution. A direct but lengthy calculation shows that S T is a growing function 
of temperature. This is in contrast to the behavior of the string tension cr T discussed in section 
3.1. Thus, the entropy contributions do play a major role at finite temperature. 

We conclude the discussion with a few short comments: 

(i) A temperature growth of the singlet internal energy at large distances was observed in the 
lattice calculations of [I]. 

(ii) The use of thermal AdS as the background metric in the low temperature phase results in 
zero entropy!^ This leads to the picture that is completely inconsistent with physics of a pure 
SU (3) gauge theory at finite temperature^ 

(hi) In Fig. 8 we plot S T versus temperature as provided by the expression (|3. 13j) . 




Figure 8: Tension E T in units of a versus temperature in units of T\. 

(iv) It is tempting to see to what extent our predictions for the low temperature behaviors of 
<j T and E T are in agreement with the lattice data of [7] . Unfortunately, the lattice data are only 

14 In this case / = 1, so (|3.10[) gives zero. 

15 Note that the Hawking-Page transition between the thermal AdS space and the Schwarzschild black hole was 
discussed in [TT] as a possible dual description of the deconfinement phase transition for large N c gauge theories. 



11 



available for temperatures close to its critical valuer I 
3.3 Polyakov Loop 

As noted earlier, at very large separation the quark and anti-quark become free. In this case the 
dominant exponent is given by Soq. So, we are in a situation in which the approximation (|2.1|) 
might be applicable. When we use it to find the expectation value of the Polyakov loop, we get 

L = e"5 5o °, (3.14) 

with the minimal area given by (|2.18p . Our normalization is stated in (|3.5p and f)3.6|) . With 
this choice, the Polyakov loop expectation value takes the form shown in Fig. 9. 




Figure 9: Polyakov loop expectation value versus temperature in units of T c . The dashed line denotes 
the value Loo. 

Certainly, it has the desired behavior: The expectation value of L is zero at low temperatures, 
while it is nonzero at high temperatures r^l One sees there is a phase transition from the confined 
phase to the deconfined phase, as expected. Numerically, the transition temperature is of order 

T c ~ 0.11\/c ~ 100 MeV . (3.15) 

This is roughly twice smaller than the value found in [T] from the spatial string tension. Thus, 
the numerical consistency is not good enough at this point. 

We chose the normalization constant in the form c(T) = aT + C. The meaning of the 
coefficients is the following: a specifies the value of L at T = oo 

Loo = e^ s - a) . (3.16) 

C specifies the form of L. Indeed, the expectation value of the loop is a continuously growing 
function of temperature for C > 0, while it has a local maximum for C < 0. The position of the 
maximum is given by a solution to equation 

CT+^ f dx exp{ = . (3.17) 

n Jo 

16 O.A. thanks O. Kaczmarek and P. Petreczky for a discussion of this issue. 
17 We consider an absolute value of L. 
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For our set of parameters a numerical analysis of (|3.17p results in T ~ 3.8 T c . 

Finally, we can compare the temperature dependence of the Polyakov loop expectation value 
as provided by our model with the lattice results of [7j- Here there is a subtle point. The 
lattice data are only available for the range 1.03 < < 6. This makes it difficult to see what 
exactly happens in the low temperature phase as well as for high temperatures. We now fix the 
normalization by fitting (|3.14p to the data given in Table 1 of [7] near T = 6T C . In Fig. 10 we 
have plotted L against We find that the temperature dependence is in good agreement for 

T>2TM 




Figure 10: Polyakov loop expectation value versus temperature. Here c(T) = 0.96 T — 0.18. The dots 
denote the data from [7j. 



4 Discussion 

What we have learned is that the 5-dimensional effective model we proposed to study the proper- 
ties of a heavy quark-antiquark pair at finite temperature turns out to be remarkably consistent 
with the qualitative expectations of thermal gauge theory. Moreover, in certain cases it provides 
the analytic results which are in good agreement with the lattice data. It is also worth not- 
ing that our model predicts the behavior of the thermodynamic quantities for low temperatures 
where field theory is unreliable and lattice data are missing. 

There are many issues that deserve to be further clarified. Let us mention some of them that 
are seemed the most important to us. 

Apparently, the model suffers from a lack of numerical self-consistency: the value of the 
critical temperature as seen from the spatial string tension [TJ turns out to be at least twice 
bigger than the value we found from the analysis of the Polyakov loop in section 3.3. There are 
two lines of thought on this problem. 

The first is to account for one- loop corrections in the world-sheet path integral (jl.5j) P^l This is 
a complicated problem. Indeed, it is a challenge to theorists to find the world-sheet formulation 
of string theory on warped geometries like AdS spaces. Among many things, it requires world- 
sheet fermions and even a world-sheet theta angle that is a two-form field B with an arbitrary 
value of J B |12j . But if it was resolved, it could help us to sum the series (jl.6p and hence refine 
the estimates. There is another interesting problem here. What were considered in [13] are some 

18 Note that a better agreement is obtained by shifting the plot of Figure 10 a little bit to the left. 
19 The model under consideration is an effective theory. It already includes some, but not all, quantum corrections 
in the approximation we used. The remaining corrections are due to string fluctuations. 
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low temperature corrections to the string tension. It will be interesting to find similar corrections 
in our model. 

The second line of thought is to somehow modify the background. It may include a slight 
revision of the metric (jl.ip or more radical changes like additional background fields. For in- 
stance, a recent proposal of [H] is even to include the tachyon background. If a modification is 
made, it would produce a number of additional free fit parameters that makes it less attractive 
for phenomenology. To escape the problem, a clever mechanism for reducing the number of 
parameters must be invented. 

Part of the interest of AdS/CFT(QCD) stems from attempts to understand the physics of 
RHIC. The reasons are the following: First, a quark-gluon plasma is strongly coupled, so the 
perturbative QCD is of limited utility. Second, the lattice data do not always provide good 
information like, for instance on transport properties. On the other hand, AdS/CFT(QCD) 
offers an alternative way of dealing with this real-world problem^ It would be particularly 
interesting to see if our model can shed some light on this subject. 

Acknowledgments 

O.A. would like to thank O. Kaczmarek, P. Petreczky, M.I. Polikarpov, and E. Shuryak for 
useful communications and conversations. The work of O.A. was supported in part by Max- 
Planck-Gesellschaft, Deutsche Forschungsgemeinschaft, and Russian Basic Research Foundation 
Grant 05-02-16486. 



References 

[1] O. Andreev and V.I. Zakharov, The Spatial String Tension, Thermal Phase Transition, and 
AdS/QCD, |hep-ph/0607026l 

[2] O. Andreev, Phys.Rev.D 73, 107901 (2006). 

[3] L.D. McLerran and B. Svetitsky, Phys. Lett.B 98, 195 (1981). 

[4] O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow, Nucl.Phys. Proc.Suppl. 129, 560 (2004). 

[5] The literature on the Wilson loops at finite temperature is vast. The following is an incomplete list: 
E. Witten, Adv.Theor.Math.Phys. 2, 505 (1998); 
S.-J. Rey, S. Theisen, and J.-T. Yee, Nucl.Phys.B 527, 171 (1998); 

A. Brandhuber, N. Itzhaki, J. Sonnenschein, and S. Yankielowicz, JHEP 9806, 001 (1998); 
D.J. Gross and H. Ooguri, Phys.Rev.D 58, 106002 (1998); 
H. Dorn and H.J. Otto, JHEP 9809, 021 (1998); 

S.A. Hartnoll and S. P. Kumar, Phys.Rev.D 74, 026001 (2006); 



S.-J. Sin and I. Zahed, Ampere's Law and Energy Loss in AdS/CFT Duality, hep-ph/0606049 
H. Boschi-Filho, N.R.F. Braga, and C.N. Ferreira, Phys.Rev.D 74, 086001 (2006); 



K. Kajantie, T. Tahkokallio, and J.-T. Yee, Thermodynamics of AdS/QCD, hep-ph/0609254 



H. Boschi-Filho and N.R.F. Braga, AdS/CFT correspondence and strong interactions, 



hep-th/0610135 



[6] D.A. Kirzhnits and AD. Linde, Phys.Lett.B 42, 471 (1972); 
L. Dolan and R. Jackiw, Phys.Rev.D 9, 3320 (1974); 
S. Weinberg, Phys.Rev.D 9, 3357 (1974). 



3 For a recent review, see [15] 



14 



O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow, Phys.Lett.B 543, 41 (2002). 
O. Andreev and V.I. Zakharov, Phys.Rev.D 74, 025023 (2006). 

E. Eichten, K. Gottfried, T. Konoshita, K.D. Lane, and T.-M. Yan, Phys.Rev.D 17, 3090 (1978); 21, 
203 (1980). 

I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, 1994. 
E. Witten, as cited in [5]; 

CP. Herzog, A Holographic Prediction of the Deconfmcmcnt Temperature, hep-th/0608151 
See the discussion of E. Witten in [5]. 

R.D. Pisarski and O. Alvarez, Phys.Rev.D 26, 3735 (1982). 

C. Csaki and M. Reece, Toward a Systematic Holographic QCD: A Braneless Approach, 
|hep-ph/0608266| 

[15] S. Gubser, Relativistic heavy ion collisions and string theory, lectures given at the PiTP 2006 program, 



http://www.admin.ias.edu/pitp/2006files/Lecture notes/Gubser Lecture notes REVISED.pdf. 



15 



